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Abstract: We theoretically and numerically investigate the capability 

of disordered media to enhance the optical path length in dielectric slabs 
and augment their light absorption efficiency due to scattering. We first 
perform a series of Monte Carlo simulations of random walks to determine 
the path length distribution in weakly to strongly (single to multiple) 
scattering, non-absorbing dielectric slabs under normally incident light 
and derive analytical expressions for the path length enhancement in these 
two limits. Quite interestingly, while multiple scattering is expected to 
produce long optical paths, we find that media containing a vanishingly 
small amount of scatterers can still provide high path length enhancements 
due to the very long trajectories sustained by total internal reflection at the 
slab interfaces. The path length distributions are then used to calculate the 
light absorption efficiency of media with varying absorption coefficients. 
We find that maximum absorption enhancement is obtained at an optimal 
scattering strength, in-between the single-scattering and the diffusive 
(strong multiple-scattering) regimes. This study can guide experimentalists 
towards more efficient and potentially low-cost solutions in photovoltaic 
technologies. 
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1. Introduction 

Maximizing the light absorption efficiency of dielectric media is paramount especially in pho¬ 
tovoltaic technologies. In general, enhanced absorption is provided by engineering the sur¬ 
face and/or the volume of the absorbing medium as to increase the time that light interacts 
with it. While attention has mainly been driven towards thin (wavelength-scale) films in re¬ 
cent years, where strong coherent effects from photonic and/or plasmonic nanostructures can 
be exploited am thicker films, such as the standard silicon photovoltaic cell, luminescent 
solar concentrators |5]|6), or dye-sensitized solar cells |7), are widely used in practice. The 
standard strategy to enhance the absorption efficiency of thick dielectric films consists in de¬ 
signing structured surfaces (e.g., gratings, random roughness, scattering particles, etc.) as to 
reduce the reflection caused by the refractive index contrast at the interface, and increase the 
optical path length in the medium by enabling an efficient light coupling to long trajectories 



in the medium IfS MTSI . see Fig. Ufa). A viable alternative to provide longer optical path lengths 
is to exploit the volume scattering provided by disordered media lfT3UT6ll . see Fig. [Q b). It is 
well known that multiple light scattering leads to longer optical paths - this property has been 
exploited, for instance, to enhance the interaction of light with gas in porous media for spec¬ 
troscopy purposes El - yet short optical paths, which dominate in reflection, tend to reduce 
the overall absorption of the medium. The respective contribution of short and long optical 
paths evidently depends on the absorption strength of the medium, suggesting that an optimal 
scattering strength, obtained by tuning the size and/or density of the scattering particles, should 
be used to reach maximum absorption enhancement. This optimum has recently been observed 
in numerical simulations on dye-sensitized solar cells ED- To date, however, the delicate and 
fundamental relation between the full path length distributions, including short and long tra¬ 
jectories, and the light absorption enhancement in weakly to strongly scattering media has not 
been investigated thoroughly. 



Fig. 1. Illustration of possible approaches to enhance light absorption in dielectric slabs, 
(a) Random roughness is introduced on the surfaces of the slab, leading to an efficient 
spreading of light in the medium. At broad angles, light follows long trajectories, (b) Vol¬ 
ume scattering is exploited to create a broad distribution of path lengths in the medium, 
consisting of trajectories that may be smaller or much larger than the sample thickness. 


In this paper, we theoretically and numerically investigate the capability of disordered media 
to enhance the optical path length in dielectric media, from the weak to the strong scattering 
regime, and maximize the light absorption efficiency. We perform a series of Monte Carlo sim¬ 
ulations of random walks to determine the path length distribution in disordered media and 
derive analytical expressions for the path length enhancement in these two limits. An inter¬ 
esting result is the fact that large path length enhancements are obtained in both weakly and 
strongly scattering limits. Because the path length distributions are very different in these two 
cases, however, the light absorption efficiency is found to be markedly different. Considering 
media with varying absorption coefficients, we find that maximum absorption enhancement is 
obtained at an optimal scattering strength, in-between the single-scattering and the diffusive 
(strong multiple-scattering) regimes. 

2. Path length enhancement 

We consider a dielectric slab, infinite in lateral directions, with effective refractive index and 
thickness L embedded in an environment with refractive index n \. The disordered medium com¬ 
posing the slab is assumed to be statistically homogeneous and isotropic, such that ballistic light 
is attenuated following the Beer-Lambert law over a typical distance that is the scattering mean 
free path £ s . For simplicity, in the analytical derivations below, single scattering is assumed to 
be isotropic, making the so-called transport mean free path £ t equal to £ s . Furthermore, light is 












assumed to be at normal incidence on the slab and completely unpolarized, making it possible 
to use average Fresnel coefficients for the reflection from the slab interfaces. 

2.1. Monte Carlo simulations 

In order to get insight into the potential of volume scattering for path length enhancement in 
dielectric slabs, we start by performing a series of Monte Carlo (MC) simulations of random 
walks in scattering media (ED- This standard numerical approach is equivalent to modeling 
light transport by the radiative transfer equation l20l . It is typically valid under the assump¬ 
tion that the scattering mean free path t s is much larger than the wavelength in the medium 
A ( kl s 1 with k = 2n/X) to avoid the occurrence of mesoscopic phenomena like localiza¬ 
tion ED, which is the case for most standard disordered media. For a sufficiently large number 
of simulated random walker trajectories, the density of random walkers at point r and time t is 
expected to be directly proportional to the light energy density averaged over disorder realiza¬ 
tions at r and t. 

The step lengths were randomly chosen from an exponential distribution, describing the 
Beer-Lambert law attenuation, with mean £ s , and the new direction of propagation after a scat¬ 
tering event was randomly and uniformly chosen on the 4 n solid angle. Reflection at the slab 
interfaces was implemented by calculating the Fresnel coefficients depending on the angle of 
incidence and averaging the reflectance over polarization. Thus, when the angle of incidence 
was greater than the critical angle, random walkers would experience total internal reflection. 

MC simulations were performed for a medium with 112 = 1.525 embedded in air (n 1 = 1) 
and for varying optical thicknesses L/l s , covering about 6 orders of magnitude (note that op¬ 
tical thickness generally refers to the quantity L/t t \ here £ t = £ s ). This system was chosen 
in view of future proof-of-principle experiments, which could be performed on polymer media 
containing dielectric (e.g. titania) nanoparticles and dyes, since ah quantities could be very well 
controlled (solutions of monodisperse titania particles are readily available and often used in 
scattering experiments, as in Ref. 11221 1. We simulated more than 10 6 random walk trajectories 
in each case and recorded the total path length I of each trajectory. Only random walkers that 
entered the medium were considered (those reflected by the first interface have a path length of 
exactly 0). Figure [2] shows the path length distributions P{1) - the integral of which equals 1 - 
in optically thick (L/£ s = 100) and thin (L/£ s = 0.0125) media, which are found to be markedly 
different. In the optically thick medium, the path length distribution is broad and smooth. This 
is a characteristic of efficient multiple light scattering, which creates trajectories with many dif¬ 
ferent path lengths, from trajectories much shorter than the sample thickness L (hence, random 
walkers necessarily escape from the illuminated interface) to very long trajectories, see the in¬ 
set of Fig. |D a). By contrast, in the optically thin medium, the path length distribution exhibits 
very sharp peaks in correspondence with multiples of L. This shows that, as expected, most ran¬ 
dom walkers escape the medium without being scattered at all. More interestingly, simulations 
indicate that a few random walkers have performed extremely long trajectories, orders of mag¬ 
nitude longer than L. This possibility is provided by total internal reflection: random walkers 
that are scattered in angles above the critical angle experience multiple internal reflections until 
they are scattered again, after an average distance £ s L, see the inset of Fig. |2jb). As we will 
show below, these rare events are those that contribute the most to the path length enhancement. 

2.2. Theory 

We now turn to the derivation of analytical expressions for the average path length in slabs 
of scattering media, starting with the case of an optically thick medium, which has been in¬ 
vestigated previously (e.g. l23l ). When the sample thickness exceeds several transport mean 
free paths (typically L > 8 £ t 1241 ), the ensemble-averaged energy density ii(r.t) can be well 



Fig. 2. Path length distributions P(l/L) (normalized to the sample thickness) in media with 
different optical thicknesses, (a) Strong multiple scattering in optically thick media (here, 
L/l s = 100) leads to a broad and smooth distribution of path lengths. Some trajectories are 
shorter than the sample thickness, and others are much longer, see the inset, (b) In optically 
thin media (here, L/£ s = 0.0125), most of the light leaves the slab without being scattered, 
but the small portion of scattered light can perform very long trajectories thanks to total 
internal reflection, see the inset. 


described by a standard diffusion equation 


du(r,t) 

dt 


= OV 2 it(r,f), 


(1) 


where D = vi t /3 is the diffusion constant, and v = cfni is the energy velocity in the medium. In 
a bounded medium, u(r,t) can be written as an eigenfunction series, the coefficients of which 
depend on boundary conditions. Considering that the randomization of the normally incident 
light takes place at a depth of one transport mean free path l t , one can solve the diffusion 
equation for a point source at a depth £ t , and reach analytical expressions for the time-resolved 
transmission and reflection 112511261 
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Here, z e = j 375 ^ is the so-called extrapolation length and r, is the average internal reflection 
coefficient, which can be calculated from Fresnel coefficients Il27l . 

The total transmission and reflection are therefore given by 


Ttnt — 


2 “ 1 . 
- L - sin 


r' co 

/ T(t)dt 

Jo 

r°° 2 " 1 

R to , = / R(t)dt = — V - sin Ml 
Jo 71 n L 


it +Z e 


L + Ze 

Ml - 

^ + 2 Ze _ 

COS 

Ml - 

L + 2 Ze _ 


(4) 


it J~Ze 
L + 2 Ze 


COS 


Ml- 


Ze 


L + 2 Ze 


(5) 




















































and the mean transmitted and reflected times of flight by 


(tr) 
{ f R ) 


Jo tT(t)dt 
Ttot 

Jo tR(t)dt 

Rtnt 


2{L + 2z e ) 2 f. 1 

) sin 
" n 3 

n= 1 


U' DTun 

2{L + 2z e ) 2 
71'DR,, 


it ~\~Ze 


T + Ze 

nn -—-— 

£ ~b 2 Ze _ 

COS 

nn --— 

£ + 2 Ze _ 


1 


■ sin 


«=1 


/Z7T- 


£t 


L + 2z r 


cos 


nn- 


Z e 


L + 2 Ze 


( 6 ) 

(7) 


The mean exit time (t) is a weighted sum of the mean transmitted and reflected times of 
flight, yielding 
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To the lowest order in £ t /L (i.e. for optically thick media) and performing the sum over n, the 
expression of (t) reduces to 
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Finally, writing ( t) explicitly in terms of £ t , a simple analytical expression for the average path 
length (I) =v(t) is found 
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In the case where the refractive indices of the scattering and embedding media match {n \ = 
n 2 ), internal reflection is suppressed (r, = 0), yielding an average path length (/) = 2.5 L. For 
a medium with refractive index n 2 = 1.525 in air («i = 1.0), we obtain r,- = 0.5841 and an 
average path length (l) — 5.31 L, which is quite significant considering the low refractive index 
contrast. Finally, it is interesting to note that the average path length is actually independent of 
l t (or is). This point will be discussed lateron. 

Let us now consider the case of an optically thin medium, L/£ t <C 1. For simplicity in the 
analytical treatment, we assume isotropic scattering, £ t = £ s . As suggested by the MC simu¬ 
lations, only a tiny portion of the incident light is scattered in the medium, yet this leads to 
extremely long trajectories in the slab due to total internal reflection, which can potentially 
result in significant path length enhancements. 

The average path length performed in the medium can be described analytically using a 
simple iterative scheme, that describes successive light scattering and outcoupling events 
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+ (1 -N Ub )P out £ 0 ut + (1 -N Ub ){ 1 ~Pout)£in 

+ (1 ~N iib )(l -Pout)Pou,£out + (1 -A^)(l Pout) 2 £in 

+ ... 

+ (1 -N iJb )( 1 -Po Ut ) n ~ l Pout£out + (1 -N ttb )(l - Pout) n £in ■ (12) 


Here, the first line describes the light that enters the system at normal incidence, the second line 
the light that has been scattered once, the third line the light that has been scattered twice, etc. 
N L jy is the portion of incident light that leaves the system without being scattered (but that can 
be reflected on the medium interfaces one or more times) and the average distance made by 































it. The portion of light 1 — No, performs an average distance s in the medium until the first 
scattering event. P ou t is the portion of light that may couple out from the slab after a scattering 
event because its propagation angle lies below the critical angle for total internal reflection. 
£ ou , is the average distance made by it, from the scattering event until it is either scattered or 
couples out from the slab. The portion of scattered light 1 — P out , instead, is guided by total 
internal reflection and performs an average distance £,„ until the next scattering event. The sum 
over n scattering events can be simplified by the power series {) a" = (1 — x) _1 , yielding 

(l) = Nijbiij, + (1 - N iib )li, s + (1 - N ub )£ out + (1 - N ijb ) ~~~~£in ■ (13) 

"out 

Let us now calculate the various quantities in Eq. (IT3l) . The attenuation of the intensity along 
the light trajectory is due, on one hand, to scattering and, on the other hand, to outcoupling 
at slab interfaces. The attenuation of the intensity (normalized to the incident intensity) due to 
scattering only follows Beer-Lambert law 

/,(/)= exp[-//4]. (14) 


Regarding the attenuation due to outcoupling at slab interfaces, we consider that light propa¬ 
gates at an angle 0 from the normal direction. At each interface, a portion of light may couple 
out, depending on the reflectance r(0), calculated from Fresnel coefficients in the case of un¬ 
polarized light. Thus, as light propagates in the waveguide, it is attenuated by a factor r along 
the path length with a period L/cos(0). To simplify the task, we assume that the decay of the 
intensity due to outcoupling may be approximated well by a simple decaying exponential 


I 0 (l) = exp 


/(I -r(0))cos(0) 
L 


(15) 


This approximation is expected to be valid to describe quantities involving many reflections on 
the interface. 

At normal incidence (0=0) and writing r(0 = 0) = /'n, the ratio of light that couples out of 
the slab without being scattered is then 
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The integrals in the denominator are present for normalization purposes (the decaying inten¬ 
sities are not probability density functions). We clearly see that for L/£ s -C 1, N,j, 1 and 
£i,b —^ (1 — ^"o)- 

Regarding the light that is scattered in the medium, the critical angle for total internal reflec¬ 
tion, 9 C , allows distinguishing two cases: 


1. 0 < 0 < 0 C : light can couple out from the slab. In the limit where £ s 2> L, the attenuation 
will be essentially due to outcoupling, i.e. scattering can be neglected in the estimation 

Of tout- 

2. 0 C < 9 < 7r/ 2: light is confined to the slab by total internal reflection. The attenuation is 
entirely due to scattering, yielding £i„ = £ s . 










The critical angle is simply defined as 
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such that the portion of scattered light that may couple out from the slab is given by 



(19) 


Here, the factor 2 comes from the fact that light can escape from both sides of the slab. 

Let us now consider the limit L/£ s <C 1 in Eq. ( 1 1 31 ). As shown above, the first term will 
converge to a non-zero value, L/{ 1 — ro). In the second and third terms, (1 — IV,/,) goes to 0, 
while both £,- s . and £ out are expected to be on the order of L. Indeed, the intensity rapidly decays 
due to outcoupling for realistic values of the reflection coefficient (an index contrast of 1.525, 
for instance, yields ro = 0.0432). Thus, in the former case, scattering of the incident light is 
more likely to occur within the first multiple reflections between interfaces, and in the latter 
case, the average distance made by reflection between interfaces for angles near normal (below 
the critical angle) is very unlikely to exceed a few sample thicknesses. As a result, the second 
and third terms can be neglected. Finally, the last term is also proportional to (1 — Ni b) but it is 
totally compensated by = £ s , which is much larger than L. 

Following these considerations, inserting the derived expressions into Eq. (Il3l) with ro = 
(, n\ — « 2 ) 2 /(ui + n 2 ) 2 , and letting l s tend towards infinity, we obtain an expression for the 
average path length 



( 20 ) 


When n\ = m, (I) = L, i.e. there is no enhancement of the path length. This is expected since 
total internal reflection is absent and cannot compensate for the tiny scattering. On the other 
hand, when n\ = 1.0 and ni = 1 .525, we obtain (l) = 4.27 L, which is surprisingly large consid¬ 
ering the vanishingly small scattering. Thus, we have shown here that, although scattering can 
be extremely weak - the slab is essentially transparent - the path length enhancement can be re¬ 
markably large thanks to the rare but very long trajectories sustained by total internal reflection 
in the slab. 

It is worthy of notice that Eq. (l20t . like Eq. ( ITTb . is totally independent of the scattering 
strength of the medium. This observation has, in fact, deep physical roots. The optical path 
length enhancement in a medium, within the radiative transfer picture, equals the light intensity 
enhancement. Upon average over incident angles, it is expected that the latter be bounded by the 
density of states ratio between the medium and its environment, i.e. by [nj/u]) 2 in 3D (yielding 
the so-called Yablonovitch limit l28l ), and be independent on microscopic details. This invari¬ 
ance property has recently been demonstrated theoretically and numerically for wave scattering 
in disordered media (29). In view of the above considerations, it is therefore not surprising that 
(l) in Eqs. ifTTb-(f20l> is independent of £ s or £ t . Because illumination is not isotropic in our 
case, different expressions are obtained in the limits of optically thin and thick media, but we 
expect that both should converge to a unique expression upon average over incident angles. In 
particular, larger incidence angles are expected to increase the average path length in optically 
thin media due to the longer trajectories provided by the unscattered light and to decrease it 
in optically thick media due to the fact that the incident light will be scattered at shallower 
depths, thereby augmenting the respective contribution of short trajectories in reflection. It will 
be interesting to investigate this aspect in a future work. 







2.3. Comparison between Monte Carlo simulations and theory 

Let us now compare the average path length obtained from MC simulations for different values 
of the optical thickness L/l s with theoretical predictions in the limits of optically thin and thick 
media. Monte Carlo simulations, performed with more than 10 6 random walk trajectories for 
each optical thickness, have been repeated 5 times in order to provide an error on the average 
path length estimation. Results are shown in Fig. [3] where the average path length is found to 
increase smoothly with the optical thickness. The larger error with decreasing optical thick¬ 
ness comes from the fact that random walkers actually experiencing a scattering event, which 
contribute for most of the path length enhancement, are fewer, making the statistics poorer. 
Nevertheless, the agreement with theoretical predictions, Eqs. CD-©. is very good, thereby 
validating our theoretical model and our understanding of the origin of large path length en¬ 
hancements in scattering media. 
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Fig. 3. Path length enhancement (/}/L determined from MC simulations (black dots) versus 
the optical thickness L/£ s , ranging from optically thin to optically thick media (left to right). 
The gray colored region indicates region where diffusion theory applies, the red and blue 
lines indicate the path length enhancement expected in the limit of optically thick (L/£ s ~S> 
1) and optically thin ( L/£ s <C 1) media. MC simulations are in very good agreement with 
theoretical predictions, Eqs. GD-Goli. 


3. Absorption enhancement 

At this stage, we can investigate the effect of volume scattering and the resulting path length 
enhancement on absorption. Absorption is assumed to be homogeneously distributed through¬ 
out the volume of the slab and is characterized by an absorption mean free path In the 
random walk picture, due to absorption, each trajectory is attenuated exponentially following 
Beer-Lambert law. The broad distribution of path lengths makes that different trajectories will 
contribute differently to the overall absorption of the medium. The medium absorption A is re¬ 
lated to the path length distribution P(Z), calculated previously and the integral of which equals 










1, as follows 



i J 


( 21 ) 


The prefactor (1 — ro) takes into account the fact that a portion of the incident light is directly 
reflected by the first interface before entering the medium. 

The absorption of the bare slab (i.e. a uniform medium containing no scatterers) can be 
derived analytically from its path length distribution P b {l) 


P b (/) = (1 - r 0 ) [8 (/ - L) + r 0 8 (/ - 2L) + ... + rgS (/ — (n + 1 )L )\, 
yielding a general expression for the total absorption enhancement T 

1 — JP(l) exp (—l/£j)dl 


( 22 ) 


(23) 


[1 — exp [—Lj £;)] / [1 — r 0 exp {—L /£;)] ’ 


In weakly absorbing media, it is known that the absorption enhancement is essentially 
bounded by the path length enhancement. Assuming that £, 2> L and £; (l), Eq. ( 12Tb reduces 

to 



(24) 


as expected. Note that the actual path length in slabs containing no scatterers, excluding light 
reflected on the first interface, equals L/(l — ro), as shown above (£,■ b tends to L/{ 1 — ro) when 
l s » L). 

Equation (f23l) was used to calculate the absorption enhancement with varying absorption 
strengths {L/ii = 1,0.1,0.01) by using the path length distribution retrieved previously in MC 
simulations on a wide range of scattering mean free paths. The enhancement factors are sum¬ 
marized in Fig.[4} a,). The first observation is that the absorption enhancement is naturally higher 
for weakly absorbing media, approaching the upper bound in Eq. (l24l) . This is expected as the 
absorption enhancement then benefits from the entire path length distribution and especially 
long trajectories, which are not immediately attenuated. More interestingly, the absorption en¬ 
hancement depends strongly on the scattering strength of the medium: (i) For optically thin 
media, T tends to remain close to 1. This is due to the fact that the rare, very long trajectories, 
typically longer than £ s , are more rapidly attenuated due to absorption (here, /j < £ s ), thereby 
not contributing fully to the absorption enhancement, (ii) In an intermediate regime, typically 
for 1 < L/i s < 10, the absorption enhancement reaches a maximum. In this optimal regime, 
scattering is strong enough to reduce significantly the amount of light that crosses the medium 
without being scattered, but not too strong to allow light penetrating the medium at a sufficient 
depth compared to L. (iii) For optically thick media, although long trajectories exist and con¬ 
tribute to the absorption enhancement, most of the light follows trajectories that are comparable 
or shorter than L, exiting the slab from the illuminated surface, thereby preventing light from 
being absorbed. When both absorption and scattering are strong, this process can dominate and 
lead to a reduction of the overall absorption, as shown for the strongest absorbing medium in 
Fig. Ha). 

The ability of light to efficiently penetrate and explore the slab volume is illustrated in 
Fig. |4jb), where the total reflection and transmission of a non-absorbing scattering medium 
are plotted as a function of the optical thickness. The optimal scattering regime clearly lies in 
the range where the same amount of light is transmitted and reflected. With increasing absorp¬ 
tion strength, L/f,-, the optimal scattering strength shifts towards lower values: light is more 
rapidly attenuated, making the presence of very long trajectories essentially unimportant. 

Finally, the actual absorption of slabs, generally more relevant than the absorption enhance¬ 
ment for applications, is reported in Table |T] for different scattering and absorption strengths. 
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Fig. 4. (a) Effect of volume scattering on light absorption in dielectric slabs for three dif¬ 
ferent absorption strengths L/ii = 1,0.1,0.01 (black, orange and blue dots, respectively). 
The gray dots indicate the upper bound for light absorption enhancement, according to 
Eq.®. The red dashed line corresponds to media containing no scatterers (enhancement 
equals 1 by definition). The absorption enhancement strongly depends on the optical thick¬ 
ness of the medium, exhibiting an optimal regime where the enhancement is maximum, 
(b) Total transmission (red points) and reflection (black points) through disordered media 
versus their optical thickness. Transmission and reflection are comparable in an intermedi¬ 
ate regime between single scattering and diffusion (strong multiple scattering). It is in this 
scattering range that light absorption enhancement is found to be maximum. 
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For the more strongly absorbing medium (L/l t = 1), the absorption enhancement is about 30%, 
leading to a total absorption of 85%, which is rather high and competitive with approaches 
relying on surface engineering |30|. 


Table 1. Total absorption in a bare slab and in slabs with varying optical thicknesses (L/t s ), 
chosen from the optimal scattering regime, and absorption strengths ( Lji ,). 



L/ti = 1 

L/li = 0.1 

II 

o 

o 

Bare slab 

0.642 

0.099 

0.010 

II 

0.794 

0.295 

0.044 

L/t s = 3 

0.847 

0.337 

0.048 

L/£ s = 5 

0.828 

0.338 

0.049 


4. Conclusions 

In summary, we have shown that volume scattering constitutes a simple, efficient and poten¬ 
tially low-cost approach to enhance the optical path length in dielectric slabs and augment their 
light absorption efficiency. We have first demonstrated via Monte-Carlo simulations and theory 
that large path length enhancements can be achieved in both optically thin and thick media. 
In the former case, this property is provided by rare but very long trajectories sustained by 
total internal reflection, contrary to the latter case, where long trajectories are due to multiple 
scattering. We have then calculated the absorption enhancement expected in scattering media 
with varying optical thicknesses and absorption strengths, and revealed the presence of an op¬ 
timal scattering regime, in between the single and strong multiple scattering regimes, where 
absorption enhancement is maximized. This optimal situation relies on an interplay between 
the actual absorption efficiency of the bare slab (characterized by L/i[) and the possibility for 
light to efficiently explore the sample volume by scattering. 

In practice, volume scattering can be obtained by embedding small scatterers (dielectric 
nanoparticles or even air bubbles) in the dielectric medium. Our study shows that an optimal 
scattering strength, driven by the scattering efficiency of the individual scatterers and the scat- 
terer density, should be used to enhance absorption efficiency, thereby guiding experimentalists 
towards more efficient and potentially low-cost solutions in photovoltaic technologies. This 
information is particularly relevant for dye-sensitized solar cells, where the size and density 
of scattering titania nanoparticles can be tuned to achieve maximum light absorption without 
affecting their primary role in separation of charges mmm. 

Finally, let us emphasize that we have restricted our study to the case of disordered media 
illuminated at normal incidence, and wherein scattering is isotropic (provided by very small 
particles) and the scatterers are non-absorbing. It will be interesting in the future to extend this 
study to media composed of resonant (metallic or dielectric) particles providing anisotropic 
scattering and to illumination at oblique incidence. While anisotropic scattering may be taken 
into account by considering that the average distance after which the direction of propagating 
light is completely randomized is the transport mean free path £ t (already taken into account 
in our diffusion model for optically thick media), oblique incidence would essentially translate 
into a modification of the unscattered light (first line in Eq. (ITU) for optically thin media and 
position of the point source in the diffusion model for optically thick media). A further step 
could be to generalize this study to multilayered structures, representing real solar cells, in 
which case Monte Carlo simulations may be more manageable than theory, mainly due to the 
difficulty to deal analytically with several boundary conditions. These advances would allow 
us to get a complete picture of the problem at hand and propose some concrete designs of 










materials operating in the visible range, where absorption enhancement is greatly enhanced, 
and that could eventually be tested experimentally. 
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